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SUMMARY 

While the vibration analysis of gear systems has been developed, a systematic approach to the reduction of 
gearbox vibration has been lacking. The technique of reducing vibration by shifting natural frequencies is proposed 
here for gearboxes and other thin-plate structures using the theories of finite elements, modal analysis, and optimiza- 
tion. A triangular shell element with 1 8 degrees of freedom is developed for structural and dynamic analysis. To 
optimize, the overall vibration energy is adopted as the objective function to be minimized at the excitation fre- 
quency by varying the design variable (element thickness) under the constraint of overall constant weight. Modal 
analysis is used to determine the sensitivity of the vibration energy as a function of the eigenvalues and eigen-vec- 
tors. The optimum design is found by the gradient projection method and a unidimensional search procedure. By 
applying the computer code to design problems for beams and plates, it was verified that the proposed method is 
effective in reducing vibration energy. The computer code is also applied to redesign the NASA Lewis gear noise rig 
test gearbox housing. As one example, only the shape of the top plate is varied, and the vibration energy levels of all 
the surfaces are reduced, yielding an overall reduction of 1/5 compared to the initial design. As a second example, 
the shapes of the top and two side plates are varied to yield an overall reduction in vibration energy of 1/30. 


INTRODUCTION 

The users of geared power transmissions desire lightweight, low-cost systems with high efficiency, large 
power capacity, longevity, and low vibration and noise. Recently, the reduction of noise caused by gears seems to be 
of growing importance, as end users perceive products with low noise and vibration to be of overall high quality. A 
project team was organized at the NASA Lewis Research Center to reduce transmission noise in helicopter cabins. 

Most vibration energy of a gear system is generated at the gear mesh. The vibration energy is then transmit- 
ted to the gearbox housing through the shafts and bearings where the energy causes housing vibration and so-called 
structure-bome noise. To estimate the vibration and noise of a gear system, the exciting source of the gear mesh and 
the dynamic behavior of the system must be considered. Once the vibration of a gearbox is determined by analysis 
or measurement (refs. 1 to 5), the sound radiated from the gearbox can be calculated by applying principles of 
acoustics (ref. 6). 

In spite of the progress made toward estimating and analyzing gearbox noise, more knowledge is necessary 
to establish methods of gearbox noise reduction. For example, an effective method of vibration reduction used for 
many products is to shift the natural frequencies of vibration to avoid resonance by adding a mass or stiffener at a 
strategic location. However, the method is not easily employed, as it is difficult to foresee the proper location for the 
mass or stiffener and then estimate the shift of the natural frequency. Recently, much progress was made in optimiza- 
tion procedures (refs. 7 to 10), and their use is becoming more common and valuable. With these ideas in mind. 
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research was conducted to develop a method of reducing transmission noise by designing the housing for low vibra- 
tion via structural optimization. The focus of this work was to design the gearbox housing for minimum vibration for 
a given level of excitation at the gear mesh. Of course, controlling the excitation source by proper design of the gear 
mesh is also important. 

A method was developed to reduce the housing vibration via the finite-element method coupled with modal 
analysis and optimization techniques. First, a self-contained and compact, yet flexible, finite element computer pro- 
gram was developed for the analysis of thin-plate and shell structures. The finite element used was a three -node 
triangular element with 18 degrees of freedom. Second, modal analysis theory was applied, and the eigenvalue prob- 
lem was solved. Third, an optimization technique was introduced to the analysis procedure. 

In this work, the overall vibration energy was selected as the objective function to be minimized by optimi- 
zation. In many cases, the dynamic response or transfer function is used as an index of vibration. For example, the 
mobility, which is the velocity caused by a unit-exciting force, could be used as the objective function. However, the 
mobility indicates the response only at the selected point. Here, the overall vibration energy was selected as the 
index of vibration and was used as the objective function for this application. The thicknesses of the elements were 
selected as the design variables, and the technique of modal analysis was used to derive vibration energy sensitivity 
with respect to the design variables. Optimum values for the element thickness were then computed by a gradient 
project method (ref. 1 1) with a unidimensional search procedure under the constraint condition of constant weight 
for the structure. 

This report describes the theory of the methods, the computer program (REVISE) based on these theories, 
and its application to some sample problems involving beams, plates, and gearbox housing. The computations are all 
based on unit harmonic forces. The frequency responses, vibration energies, and optimum shapes are calculated and 
illustrated. Finally, results and conclusions from the sample problems are discussed concerning the newly developed 
noise reduction method and its application to gearbox noise reduction. 


SYMBOLS 


a 

shape function coefficient 

A 

element area 

b 

shape function coefficient 

c 

shape function coefficient 

[C] 

damping matrix 

{d} 

steepest descent vector 

D 

flexural rigidity 

E 

Young’s elastic modulus 

f 

frequency 

(fit)} 

exciting force vector 

{F} 

modal exciting force vector 

g 

behavior constraint 

h 

side constraint 

[I] 

unit matrix 

k 

stiffness 

[K] 

stiffness matrix 

L 

shape function 

(m,n) 

mode shape 

[M] 

mass matrix 

N 

shape function 

N 

degrees of freedom 

P 

load 

[P] 

projection matrix 

[Q] 

coefficient matrix 
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s 

step size 

t 

thickness 

T 

vibration energy 

[T] 

coordinate transformation matrix 

u 

displacement; displacement in x direction 

(U(t)} 

displacement vector 

U(t)} 

velocity vector 

{ ii(t)} 

acceleration vector 

U 

work done per unit load 

u* 

work done by element unit strain 

{U} 

complex amplitude vector 

V 

displacement in y direction 

w 

displacement in z direction 

w R 

rigid body motion displacement 

* 

w 

element strain displacement 

w 

weight of structure 

w 

prescribed weight of structure 

{X} 

design variable 

{y} 

undamped eigenvector 

X, Y, Z 

Cartesian coordinates 

{z} 

damped eigenvector 

a 

proportional damping factor for mass 

P 

proportional damping factor for stiffness 

§rs 

Kronecker’s delta 

{§} 

nodal displacement vector 

{5*} 

relative nodal displacement vector 

0 

element rotation 

X 

undamped eigenvalue 

[A.] 

diagonal matrix of eigenvalues 

[A] 

coordinate transformation matrix 

P 

damped eigenvalue 

V 

Poisson’s ratio 

!q(t)} 

modal coordinate displacement vector 


modal coordinate velocity vector 

P 

mass density 

w 

normalized eigenvector 

m 

modal matrix 

CO 

undamped natural angular frequency 

Subscripts 

B 

submatrix 

ex 

excitation 

n 

number of modes considered 

0 

initial 

q 

number of constraint surfaces 

r 

counter for order of eigenvector 
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rl imaginary part of r-th component 

i'R real part of r-th component 

s counter for order of eigenvector 

Superscripts 

e element 

k k-th iteration 

T transpose of matrix or vector 


THEORETICAL FUNDAMENTALS 
Modal Analysis of Forced Vibration 

Vibration analysis of linear systems based on the mode shapes of natural vibration is commonly called 
modal analysis. The main purpose for applying modal analysis here is to convert a set of simultaneous differential 
equations into a set of independent differential equations by linear transformation of coordinates. Let an elastic 
structure be discretized into finite elements with A degrees of freedom. The equation of motion can then be given in 
matrix form as 


[M]{ii(t)} + [C]{u(t)} + [K]{u(t)} = {f (t)} (1) 

where [M] is the mass matrix, [C] is the damping matrix, [K] is the stiffness matrix, { u(t) } is the displacement vec- 
tor, and { f(t) } is the exciting force vector. By assuming proportional damping, the damping matrix may be repre- 
sented by a linear combination of the mass matrix and stiffness matrix as 

[C] = oc[M] + p[K] (2) 

The eigenvalue problem that is related to equation (1) for the case of proportional damping can be de- 
scribed as 


p 2 [M]{z} + p(a[M] + p[K]){z} + [K]{z} = {0} 
The eigenvalue |1 and the eigenvector { z } in equation (3) are obtained from 


X 


(p 2 + exp) 

(Pp + 1) 


{z} = {y} 


(3) 


(4) 


where X and { y } are the eigenvalue and eigenvector of the undamped vibration system with the same matrices [M] 
and [K] such that 


[K]{y} = X[M]{y} 


(5) 


If matrices [M] and [K] are positive definite, as is usual for vibration systems, then X, { y } , p, and { z } are 
real and the eigenvector { z } has the orthogonal properties 
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( 6 ) 


{z r } T [ M ]{z r } = m r, 

{ Z r} T [ M ]{ Z s} = 0 

(r*s) 

{z r } T [ K ]{z r } = k r , 

{z r } T [ K ]{Zs} = ° 

( r *s) 

{z r } T [ c ]{ z r } = C r =am r +P k r, 

{ z r} T [ C ]{ Z s} = 0 

(r*s) 


In equation (6), subscripts r and s indicate the order of the eigenvectors. By ignoring modes of vibration 
beyond the frequency range of interest, we can reduce the degrees of freedom to be solved (i.e., r, s = 1 ... n 
(n « AO), where n eigenvectors are considered in the analysis. The related eigenvalues are given by 


A 


r 




(7) 


where co r is the r-th undamped natural angular frequency. By applying equation (6), the set of N simultaneous equa- 
tions that describes the motions of the system (eq. (1)) can be transformed and reduced into n independent 
differential equations. 

In the particular case of harmonic vibration, the transformation mentioned above can be done more simply. 
Let { tf) } be the eigenvector which is normalized with respect to mass matrix, 

MViKHrs (8) 

where 8 rs is the Kronecker’s delta. The modal matrix, [<t>], is defined by the assembly of these eigenvectors as 


[°]=[M- W.-.R}, 


(9) 


From equations (8) and (9), and applying the principle of Rayleigh's quotient, the following relation is 
derived instead of equation (6) 


[<&]>][<&] = [I] 

W T [KW = W 


( 10 ) 


In equation (10), [I] is an nxn unit matrix and [A.] denotes the diagonal matrix composed of n eigenvalues. 
By the orthogonal properties of the modal coordinates, any displacement { u(t) } can be approximately represented by 
a linear combination of the eigenvectors 

{u(t)} = EWt){*rH'l>E(t)} (11) 

r=l 

When the structure is excited by a harmonic force of angular frequency to, then the force can be described 
by 

{f(t)} = {F} e j Q)t (12) 

and the steady-state vibration of the structure can be written in the same form as 

{u(t)} = {U}e j “ t = [4>]{y e j a)t (13) 
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Note that here, { U } is independent of time. Substituting equations (12) and (13) in equation (1) and pie- 
multiplying by [<t>] gives 

—CO 2 [<&] T [M][<&]{^}e j mt + jco[0] T [C][<F]{^}e ja)t + [<I>] T [K][<D]{^}e jcot = [<F] T {F}e ja,t (14) 

By applying equation (10) to equation (14), n independent algebraic equations are obtained for the evalua- 
tion of the vector {^} in matrix form as 


-to 2 {^} + jco[a[I] + P[?i]]{^}+ [*]{£} = [<f>] T {F} 


(15) 


The r-th component is given by 




RjV) 

A, r -co 2 ) + jco(a + |3?i r ) 


(16) 


where r = 1 , . . . , n. Consequently, the complex amplitude { U } of the forced-vibration problem is obtained by substi- 
tuting {^} into equation (13). 


Vibration Energy and its Sensitivity 

The velocity { u(t)} is generally represented in the same way as the displacement (eq.(Tl)) as 

{u(t)} = [<D]{4(t)} (17) 

In case of harmonic excitation, from equation (17) the velocity can be written in the same form as 
equation (13), 

{(i(t)} = {u} e j“ t = [<l>]{4}ej a,t (18) 

The velocity can be directly derived by differentiating the displacement equation (13), directly giving 

{(i(t)} = jco[4>]{yej fflt (19) 

From equations (18) and (19), the vector { £ } is given by 

{4} = jco{^} (20) 

Therefore, the real and imaginary parts of the r-th component is obtained as follows 
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4rR “ ' 


(oc + (3A. r ){(|) r } T {F} 
V(^r - ® 2 ) +(oc + (3A, r ) 2 


— (A, r - co 2 ){(]) r } T {F} 

e _ CO v ' 

Sri “ -j 2 _ 

(A. r — co 2 ) +(a + (3A r ) 2 
_co“ v ’ 


( 21 ) 


The vibration energy T of the structure is given by 

T KjVlM 

2 

{r} T [*] T [Mj[*]{4} 

2 

(t-m 

2 

where * indicates the conjugate complex; therefore, the energy is represented by 

T = aL& + ^i) 

r=l 


(22) 


(23) 


The sensitivity of the vibration energy with respect to the design variable x; (i = 1, . . m) is obtained from 


3T 

3x i 



r=l V 


3k] 

3 Xi ^ rI 9xJ 


(24) 


The partial derivative terms of equation (24) can be determined from differentiating equation (21) by the 
design variables. Since a, [3, CO, and { F} are independent of the design variables, the sensitivity of the vibration en- 
ergy can be found once the functions that describe the sensitivities of the eigenvalues and the eigenvectors are 
known. An explanation of their derivation follows. 


Sensitivities of Eigenvalue and Eigenvector 

The derivation of the sensitivities of the eigenvalues and eigenvectors, summarized briefly in this section, is 
based on the method proposed by Fox (ref. 12). 
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The r-th (r = 1, n) order eigenvalue X,- and eigenvector { cf) r } are obtained from the following equation: 


[K]K} = MM]k} 

Differentiating equation (25) by the design variable yields 


( 25 ) 


dX i 

dx- 


[M]K} = 


3K 

VL 3x iJ 


-X, 


m 

dx- 


i J J 


{^r} + ([K]-^[M]} 


3<t> r 

dx- 


(26) 


Premultiplying equation (26) by { 4> r } T { ?e r } T , the second term on the right-hand side vanishes because 

k} X ([K] - MM]) = [([K] - MM])k}f = 0 (27) 


Therefore, the sensitivity of the eigenvalue is given by 


dX r 

3x i 


M 1 


VL 


3K 

dx- 


-x. 


3M 

3x- 


M 


jy 


(28) 


Although equation (26) includes the sensitivity of the eigenvector, it can not be solved as written because 
the determinant of the sensitivity’s coefficient matrix is equal to zero. So, an equation is introduced independently by 
differentiating the first relation of equation (10) to give 


k}>] 


dx { 




(29) 


An appropriate row of matrix equation (26) is replaced by equation (29) to avoid a zero determinant. Then, 
the sensitivity of the eigenvector is calculated by solving the simultaneous equations. The row that relates to the 
largest eigenvector should be replaced to obtain the best solution. 


Search for the Optimum Value 


The optimum design problem for this analysis is expressed as follows: 


Minimize 

Behavior constraint 
Side constraint 


Design variables 


T{x} 

g{x} = W - W{x} = 0 
h 1 {x} = x i -x min >0, 

M x } = x max- x i 

{ x } = ( X 1’ X 2 ’"-’ x n) T 


(30) 


The elements’ thicknesses are the design variables { x } . The prescribed weight and the weight of the struc- 
ture are represented by \y and W{x}, respectively. By the behavior constraint, the designs are of constant total 
weight. Minimum and maximum limits are set for the elements’ thickness by the side constraints, hi and h?. 

Many procedures search for the optimum value for this kind of problem by first obtaining the feasible di- 
rection along which the design variables are modified to minimize the objective function. In this analysis, the gradi- 
ent projection method (ref. 1 1) (which can be applied to problems with linear or nonlinear constraints) is used to find 



this direction. Following immediately is an outline of the procedure for the case of linear constraints as it applies to 
this research. 

In the feasible region, the point corresponding to the current design is moved in the direction of the steepest 
descent vector {d (k) }, and a new design is obtained by 


where 


{ x ^l)} = { x (k)} + s { d W} 


K1=- 

3T ({x <b> }) 

dx 




(31) 


In this equation, the superscript (k) means the k-th iteration, and s is the step determined by a unidimen- 
sional search. When the new point { x c } = { x (k+l '} reaches the intersection formed by the constraint surfaces, the 
steepest descent vector is projected to the intersection to obtain the next feasible direction. The projection of the 
descent vector is given by the following relations: 


where 


{d P } = [P]{d{x c | 


[P] = [I]-{U}[V]{U} T 
[V] = [{U} T {U}] 1 


{uMkl’kl’-’K} 


w= 


3xj 3 xt 3x n j 


(32) 


To simplify, the constraint functions g({x}), hi({x}), and h 2 ({x}) are written as g; in equation (32). The 
matrix [P] is called the projection matrix. In a problem with linear constraints, the constraint surfaces are all planes. 
Therefore, the point obtained from equation (31) by using the projected vector { d p } always stays in the feasible 
region. 

The steepest descent vector {d (k) } is formed by the sensitivity or the first derivative of the objective func- 
tion as shown in equation (31). If the second derivative was obtained, the step in equation (31) could be evaluated by 
a faster method (e.g., the Newton-Raphson method). However, this would require the evaluation of the second de- 
rivatives of the eigenvalues and eigenvectors, which is fairly complicated. Therefore, the step is evaluated by the 
simple procedure explained below. Since the design variables in the current design { x fk| } are known, the value of the 
objective function can be evaluated. Let Tj be the value, where the subscript 1 means the current design. Similarly, 
let To*- 1 * and T 3 ‘ 1 1 be the values of the objective function at the different states of design variables of the 
i-th iteration. The procedure can be itemized as follows: 

(1) Estimate T \ 

(2) SetT 2 (1) = T 1 ,T 3 (1, = T 1 

(3) Assume the step s klk 

(4) Evaluate T^ lk by using s (1) and equations (31) and (32), if necessary 

(5) Compare T s (1) to T 2 (1) 

(6) IfT s (1) >T 2 (1) 

then set T 3 (2) = T s (1 \ = s kl - / 2 

evaluate T s ® for T 2 (2) 

return to (5) and continue until T s (l ) < T 2 ' l ~ l) 
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(7) IfT s (1) <T 2 a) 

then set T 2 <2) = s^ 2) = 2xs^ 

evaluate T s (2) for T 2 (2) 

return to (5) and continue until T s (l) > T 2 (b 

This procedure finds the closed region [{x®}, {x (k) } + s (1) {d p }] or [{x®}, {x (k) } + 2s (b {d p }], which con- 
tains the relative minimum. The lower and higher limits of the region give the value T ] and T/' 1 , and the middle 
point give the value T 2 fl) . If the region is found, the objective function is approximated locally by a quadratic ex- 
pression based on these values, and the step that gives the minimum point of the quadratic expression can be deter- 
mined. The procedure is ineffective for finding the region when the objective function decreases monotonously. In 
such a case, the procedure is stopped after the appropriate iterations and the step at the last iteration is taken as the 
solution. The value of the energy in the last iteration is stored and used for the judgment of convergence. The judg- 
ing criteria for convergence is l/1000th relative variation. 


OPTIMUM PROGRAM DESIGN 
Finite-Element Formulation 

The best way to model the free surface of three-dimensional structures (e.g., a gearbox) for use with an 
optimization procedure may be to use triangular shell elements. Many kinds of elements have been proposed and 
successfully applied to both static and dynamic analyses. If the total degrees of freedom is constant, using a small 
number of high-precision elements generally leads to good results. However, a larger number of nodes is needed to 
describe the three-dimensional geometry of gearbox housing. Consequently, since the total degrees of freedom for 
the analysis must be practical, a three-node triangular element of 18 degrees of freedom was selected. The element 
matrices of the shell element were formulated by combining the matrices of plate bending and plane analysis prob- 
lems (ref. 13). The formulation is summarized as follows: 

Finite Element for Plate Bending . — The actual displacement of middle surface w (fig. 1) is represented by 

. . . . . 

the sum of the deflection of the supported element w , which induces the strain in the element, and the displacement, 

D 

w , caused by rigid body motion, 

w = w* + w R (33) 


Using area coordinates, where i = 1, 2, 3, 

_ (tq+biX + Cjy) 

1 A 


a i 

= x i+ ty i+ 2 -x i+2 y i+ i 


b ! 

= yi+i -Ni+2 

(34) 

C . 

= x i + 2~ x i + l 


A 

— cij + 3-2 “1” 3^ 



The displacement caused by rigid body motion is expressed as 


3 

w R =^F i w i (35) 

i=l 

where Wj denotes the displacement at node i. In equation (34), if the subscript i+1 or i+2 is larger than 3, it represents 
the value subtracted by 3. The relative displacement at node i is defined by the rotations about the x and y axes. 




3w x 3w x 
3y dx 


0 x- 


dw 


dy 'V 


dw 


R V 


dx 


(36) 


to 



The relation between the relative nodal displacement {8*} e and the actual nodal displacement { 8 } e is given 


by 


{«*}'= M8>' 


(37) 


The matrix [T] is slightly different from the one given in reference 13 because of a difference in the defini- 
tion of the nodal displacement in equation (36). The deflection w is represented by 


w*=[n 1 ,n 2 ,n 3 ](5;,5^5;) 1 

=t( N K +N K, 


i=l 


(38) 


and the relative shape functions Ny and N; x are defined as follows: 


2,1 


Nf = b- , , L; , ~L- H L;L- . ,L; , o — b : , T L- . ,Lt H — L-L- . ,L- 


2,1 


"i - u i+l| i+2 i T 2 i i+1 i+2 I u i+2l ^i+l^i ^ 0 ^i^i+l^i+2 


N i _C i+l( L i+2 L i + 9 L i L i+l L i+2l C i+2 L i+l L i + 9 L i L i+l L i+2 


(39) 


The element stiffness matrix is evaluated by 


K B ] e =m T [KBr[T] 


sfc t t 

The (i,j) submatrix of [K 3 ] is given by 

K B i ,HM k ‘] dxd >' 


k * 

11 = 


ki2 = 


k * 

21 - 


k 2 2 - 


d 2 Nf 

d 2 Nj 

dxx 

dxx 

d 2 Nf 

d 2 Nj 

dxx 

dxx 

d 2 Nf 

d 2 Nj 

dxx 

dxx 

d 2 Nf 

d 2 Nj 

dxx 

dxx 


+ v 


d 2 Nf 3 2 N x (3 2 N x 3 2 N- 3 2 Nf 5X1 3 2 Nf 5 2 N x 3 2 N x 3 2 N| 

1 “ l \ t 1 J i i J _i 1 L i on 1 J 


3xx 3yy 3yy 3xx 


^2xtX 


3yy dyy 


+ 2(1 — v)- 


dxy dxy 


+ v 


a 2 Nf 


9 N J d 2 Nf 
- + 1 


3 2 N^ a 2 Nf 3 2 N? 


dxx dyy dyy dxx 

\ 

f a2 XT y ;im x a\i x ) 


+ V 


+ v 


dyy dyy 


J on i d z N x 
— + 2(1- v) 


3 N I 


dxy dxy 


d 2 Nf d 2 N x d 2 Nf d 2 N x 

+ — — 5 — i — — + 2(1 — v) - 1 J 


dxx dyy dyy dxx 


d 2 Nf a 2 Nj d 2 Nf 3 2 Nj 
dxx dyy dyy dxx 


dyy dyy 


dxy dxy 


d 2 Nf 3 2 N[ d 2 Nf 3 2 N[ 

+ — — - — - — — + 2(l-v)- 1 J 


dyy dyy 


dxy dxy 


(40) 


(41) 
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Here, D denotes the flexural rigidity of the plate. Every term in equation (41) is calculated by substituting 
equation (34) into equation (39) and then differentiating by x and y. For example, d 2 N)/( : ) xx is given by 

3 N -=-V( b i+l-bi+2) b i+l b i+2 L i 


3xx A 


+ J2 {- 2b i b i+2 +( b i + l - b i + 2) b i + 2 b i}L i+1 (42) 

+ -^2"{ 2b r b i+l +(b i+1 -bi+2) b i b i+i}Li+2 

The integration in equation (41) includes the combination of the integration of LjLj, where i,j = 1, 2, 3, and 
has the form 

| jLjLjdx dy = + ^ b i b j(*I +4+4) 

+ ^( b i c j +b j c i)( x iyi +x 2y 2 +x 3y 3 )+^; c i c j(y? + yi + y|)} ( 43 ) 


The displacement w is rewritten in matrix form as 

(N){S}' 


w = 


(N) = (l„ Nf, Nj, L 2 , N|, N 2 , L 3 , N 3 , Nl) 

{5} = ( Wl , e: b e; b w 2 , 0* x2 , e; 2 , w 3 , e; 3 , e; 3 ) T 

The relationship between { 8 } e and { 8} e is written in the same form as equation (37): 


(44) 


8 = T 


r]{8 }' 


(45) 


The matrix [ -p] of equation (45) has a dimension of 9x9. Substituting equation (45) into equation (44), we 


obtain 


w = (N){8} e 

(N) = (n)[t] = (L)[Q] 

(L) = (Lj, L 2 , L 3 , Lj L 2 2 , L 2 L 2 3 , L 3 L 2 , L 2 L 2 , L 2 2 L 3 , L 2 3 L h Lj L 2 L 3 ) 


(46) 


where matrix [Q] is the coefficient of (L) and has dimensions of 10x9. The consistent mass matrix [Mb] 6 is finally 
evaluated by the method of numeric integration 


[M B f pt(N) T (N)dx dy = pt[Q] T [J J (L) T (L)dx dy][Q] 


( 47 ) 


where p and t are the mass density and element thickness, respectively. 
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by 


Finite Element for Plane Analysis Problems . — The shape function for the plane analysis problem is given 


N i = 


( a i + bjX + Ci y) 


The consistent mass matrix [M p ] e is evaluated by the equation 


[M p ] C =JJpt(N) T (N)dxdy 
[N] = (N 1 [I],N 2 [I],N 3 [I]) 


(48) 


(49) 


where [I] is a 2x2 unit matrix. The (i,j) submatrix of the stiffness matrix [K p ] e is written 





b lbj 


+ 


1 — v 
2 


u 1_V U 

vc i b j + ^ b i c j 


, 1_V u 

vb i c j+^ c i b j 
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The Shell Finite Element . — The nodal displacement is defined by the following expression, where i=l, 2, 3 


{8j} = (u, v, w 0 X , 0 y ,0^ 


(51) 


Here, u, v, and w are displacements along x, y, and z axes while 0 X , 0 y , and 0 Z are the rotations around x, y, 
and z axes. 

The stiffness matrix for the shell element is obtained by combining the stiffness matrices of the plate-bend- 
ing and plane analysis problems; therefore, the (i,j) submatrix has the form 
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In equation (52), {Kg z } indicates the stiffness for rotation about the axis normal to the element. However, it 
is difficult to formulate theoretically. The stiffness term {Kg z } is given approximately by the following matrix as was 
used in the present work: 
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where a = 0.03. 


( 53 ) 
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The mass matrix of the shell element can be formulated in the same way as the stiffness matrix. The trans- 
formation of the element stiffness matrix to the global coordinate system is defined in the ordinary way, 

[ K ]global=[ A ] T [ K ]local[ A ] (54) 

where the 18x18 matrix [A] is formed by putting the 3x3 direction cosine matrix in the diagonal submatrix. 


Structure of Computer Program REVISE 

The computer code REVISE (REduction of Vibration in Structural Engineering) was developed following 
the procedures described in the previous sections and its structure is shown in figure 2. The program consists of two 
main parts. The first part is the analyzer, which is composed of the finite-element analysis, the eigenvalue analysis, 
and the evaluation of the vibration energy. The second part is the optimizer, which includes the sensitivity analysis 
and the search procedure for the optimum value. Two IMSL (ref. 14) routines are used to calculate simultaneous 
equations and invert matrices. The user can control the flow of the computations by input parameters and may obtain 
the result of the eigenvalue analysis only. This control is useful to determine the dynamic behavior of the model 
before executing the optimum design. 

The program was developed on a Digital Equipment VAX cluster at NASA Lewis and was executed inter- 
actively in case of a smaller model (e.g., a beam or a plate). On the other hand, in case of a larger model such as a 
gearbox, the program was submitted to a CRAY XMP/28 supercomputer. The program is coded in FORTRAN77, 
and no particular hardware or software is required; therefore, users can modify it to introduce other methods of 
optimization, different finite elements, and so forth. 


Accuracy of Numerical Solutions 

Validating REVISE has two aspects: verification of the numerical solution method and the estimation of 
errors. Estimation of the errors requires experimental data, and some comparisons to measured values are presented 
later with the results of the gearbox optimization. However, verification of the numerical method can be accom- 
plished by comparing numerical solutions of problems to exact analytical solutions. Following are results from a few 
simple examples solved to verify the accuracy of the finite element and modal analysis solutions calculated by the 
computer program. The material constants used were modulus of elasticity E = 206 GPa, Poisson’s ratio v = 0.3, 

o 1 

mass density p = 8000 kg/m , and coefficients for defining the proportional damping matrix a = 1 .0 s and 
P = 5.0 x10“ 7 s. 

Static Bending of a Simply Supported Square Pla te. — Table I shows the numeric solutions and the exact 
analytical solutions for the maximum deflections of a simply supported square plate for two loading conditions. One 
condition is the deflection caused by a central-concentrated load P and the other is the one caused by a uniformly 
distributed load pL . The solutions for deflections converge well toward the exact values, and the error of the nu- 
meric solution is less than 1 percent even for the case a coarse mesh of 8x8x2 elements. 

Static Bending of a Cantilever Plate . — Jaramillo (ref. 15) obtained the exact analytical solutions for the 
deflection and bending moment of a cantilever plate with infinite width and a concentrated load. His result shows 
that the deflection diminishes rapidly as the distance from the loading point increases, and the deflection becomes a 
negligibly small value at a position along the width of more than three times the length. Therefore, a plate with a 
width six times its length may be considered approximately as a plate of infinite width. Such a plate was divided into 
an 8x32x2 mesh, and the numerical solutions were obtained. Comparing the numerical solution to the analytical 
solution of Jaramillo, the maximum deflections were within 1.8 percent of each other. The maximum bending mo- 
ments were within 4.7 percent of each other. 

Static Analysis of a Cylindrical Shell . — Scordelis and Lo (ref. 16) obtained the analytical solution for the 
deformation of a cylindrical shell caused by its weight. The numerical solutions are compared with their analytical 
solutions in figure 3. The maximum deflection is close to the analytical solution. On the other hand, the circum- 
ferential bending moment is about 10 percent less than the analytical solution. Such errors for the bending moment 
is typical for finite-element solutions by the displacement method as used here. 
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Eigenvalue Analysis of a Simply Supported Rectangular Plate . — A rectangular plate, 600- by 400-mm and 
5-mm thick, was divided into meshes of 4x4x2, 6x6x2, and 8x8x2 dimensions. A similar mesh pattern was later 
used to optimize a gearbox housing. The computed natural frequencies of the plate are shown in table II. In the table, 
the mode shape indication of (m,n) means the mode shape has m-1 and n-1 nodal lines along the two sides. The re- 
sults of the numerical computations are mostly smaller than the analytical solutions, with the error depending on the 
mode shape and order. To provide another comparison, the arithmetic mean of the squared error was calculated, and 
its square root is illustrated in figure 4. Based on this arithmetic mean error, if a 6 percent precision is desired, then a 
6x6x2 mesh of elements is sufficient for solutions of modes up to 10th order. 

The numerical solutions of the previously described sample problems showed good agreement with analyti- 
cal solutions, thereby demonstrating and verifying the finite element and modal analysis done via REVISE. 


OPTIMUM DESIGN OF BEAMS AND PLATES 

To illustrate the optimization process and verify REVISE, the optimum shapes of beams and plates were 
computed. The material constants used were modulus of elasticity E = 206 GPa, Poisson’s ratio v = 0.3, mass den- 
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sity p = 8000 kg/m , and coefficients for defining the proportional damping matrix a = 1 .0 s and 
P = 5.0x10“ 7 s. 


Optimum Shapes of Beams 

Because the element used in this research is not a beam but a shell element, the computer program is not 
specifically suited for beam analysis. However, a sufficiently narrow-width plate can approximate a beam, and is 
referred to as a beam in this section. The pattern of the finite-element mesh is shown in figure 5. The beam is first 
divided into a 10x2 mesh of equal rectangles, and then each rectangle is subdivided into two triangles such that the 
pattern is symmetrical. The element numbering is shown in figure 5. The beam dimensions are length L = 200 mm, 
width W =10 mm, and initial thickness to = 5 mm. 

Simply Supported Beam . — Comparing the numerical solution of the eigenvalue problem to analytical solu- 
tions, the first four natural frequencies and their errors are 288.9 Hz (0.5 percent), 1 172 Hz (1.9 percent), 2704 Hz 
(4.5 percent), and 4982 Hz (8.3 percent), respectively. With the center of the beam excited by a unit harmonic force, 
the frequency response of the vibration energy was evaluated numerically at every 10 Hz using the first four modes. 
This result is illustrated by the broken line in figure 6 with the abscissa representing the excitation frequency. The 
peak in the figure shows the resonant response of the fundamental natural mode whose frequency is 288.9 Hz. The 
beam’s vibration energy of about 0.207xl0 -7 J in response to a unit harmonic force of 270 Hz applied at the center 
is considerably high because of the influence of resonance. 

The optimum design was calculated to minimize the vibration energy when excited by the 270 Hz 
frequency force and with the upper and lower limits for the beam thickness constrained to be 1 < t < 10 mm. At each 
iteration of the optimization procedure, the beam’s response energy was calculated and normalized to that of the 
initial beam of constant thickness. The trend of the normalized response energy with iteration number is shown in 
figure 7. The energy decreases rapidly and converges to about one tenth of the initial energy after four iterations. The 
frequency response of the optimally designed beam is indicated by the solid line in figure 6. The fundamental natural 
frequency was shifted from 288.9 to 319.7 Hz. 

The optimum shape of the beam is illustrated in figure 8. The thickness is greatest at the center and gradu- 
ally decreases towards the supported ends. The shape of the optimized beam is compared to the shape of a beam of 
uniform strength in figure 9. The latter is well known as the optimum shape for maximum stiffness under static load 
(ref. 17). Though the thicknesses near the loading points differ, the overall shapes of the beams are very similar. 
Another optimum shape, obtained by using a more dense finite element mesh of 20x2x2, is also shown in the figure. 
It is approximately the same shape as obtained for the coarse 10x10x2 mesh. 

In the above optimizations, the four elements across the beam width, (e.g., elements 1 to 4 in figure 5) were 
forced to equal thicknesses. Eliminating that constraint from the optimization procedure, we obtain the optimum 
shape illustrated in figure 10. At most locations, the beam thickness is at the higher or lower constraint limits. How- 
ever, the tendency that the thicknesses decrease toward the supporting ends is the same as the former examples. The 
variation of element thickness as a function of the number of optimizing iterations is summarized in table III. In the 


15 



first few iterations, the element thicknesses across the beam width are very close to each other, then, as more itera- 
tions are completed, the thicknesses gradually expand. 

The frequency response and the convergence of the total vibration energy are shown in figures 1 1 and 12, 
respectively. As compared with the former examples of constant thickness across the beam width, more iterations are 
needed for convergence to the optimum shape. Also, the total vibration energy is reduced by a greater amount. The 
beneficial effect of the optimized beam shape is illustrated in the continuative vibration amplitude plot responses to a 
unit harmonic force, figure 13. These responses were evaluated numerically by the method described in the Modal 
Analysis of Forced Vibration section. The response amplitude was reduced by an additional order of magnitude for 
the second optimization, allowing thickness variation across the beam, compared to the first optimization with con- 
stant thickness across the beam’s width. 

As another example, the vibration energy at 2680 Hz was minimized for the same simply supported beam. 
This excitation frequency is a little lower than the third natural frequency. First, the thickness across the beam was 
constrained to be uniform and the optimized shape was calculated. Figure 14 shows the frequency responses of the 
total vibration energy before and after the optimization, and the convergence of the energy is shown in figure 15. The 
optimum shape is illustrated in figure 16. The thicknesses near the nodes of vibration are decreased by the optimiza- 
tion method. Next, the beam thickness across the width was allowed to vary, and the optimized shape recalculated. 
The thickness distribution of the optimized shape, shown in figure 17, is almost same as the first optimization, figure 
16. The frequency response and the convergence of the energy for the second optimization are shown in figures 18 
and 19, respectively. The resonance peak shifts to a much higher frequency for the second optimized shape, figure 
18, compared to the first, figure 14. The continuative vibrations of the original beam and the two optimized shapes, 
figure 20, dramatically illustrates the benefit of the optimization procedure to reduce vibrations. 

Clamped Beam . — As another sample problem, a beam with the same geometry as discussed in the Simply 
Supported Beam section, but with clamped rather than simply supported ends was studied and optimized. The first 
four natural frequencies of the initial beam with constant thickness were calculated numerically as 671.2, 1904, 

3887, and 6787 Hz. These numeric solutions are not quite as accurate as those for the simply supported beams. The 
frequency response of the initial-thickness beam was evaluated using only the first four modes. The frequency re- 
sponse is shown in figure 21. Next, the beam was optimized to minimize the vibration energy at 650 Hz, which is 
slightly lower than the first fundamental natural frequency. 

The thickness across the width of the beam was constrained to be equal for the first optimization. The con- 
vergence of the energy is shown in figure 22. More iterations are required for convergence as compared to the optim- 
ization of the simply supported beam. The optimum shape of the beam is illustrated in figure 23. The thickness 
increases slightly at the center and greater at the fixed ends. These positions correspond to positions of maximum 
bending moments. The optimum shape of the clamped beam was compared with the shape of a beam of uniform 
strength, figure 24 and the shapes were found to be similar. However, the shape optimized for minimum vibration 
has a greater thickness near the fixed ends while the uniform strength beam is thicker near the center. The difference 
seems to be effective in reducing vibration. 

A second optimization of the clamped beam was done, eliminating the constraint of equal thickness across 
the width. The optimum shape is shown in figure 25. The mean thickness for each group of four elements across the 
beam width are close to the thickness distribution of the first optimized shape. The frequency response of the vibra- 
tion energy and the convergence of the optimization are shown in figures 26 and 27. Comparing figure 21 to figure 
26, the resonance peak is shifted much higher in the case of allowing rather than constraining thickness variation 
across the beam width, leading to a lower energy level. Figure 28 shows the amplitudes of continuative vibration of 
the original and two optimized clamped beams. The amplitude is reduced greatly by the optimization. In figure 28, 
the slope at the clamped end is not zero because the amplitudes at the nodes are connected by straight lines rather 
than curved. Therefore, the slope not equal to zero is not indicative of poor solutions. 


Optimum Shape of Plates 

As another sample problem, the optimum shape of a rectangular plate was calculated. The plate size, fig- 
ure 29, was 600- by 400-mm with an initial thickness to of 5 mm. To define the finite-element mesh, the plate was 
first divided into a 6x6 set of equal rectangles. Then, each rectangle was subdivided into two triangles such that the 
pattern was symmetrical (fig. 29). The element thicknesses were constrained to be 2 < t < 10 mm for optimization. 
Before executing the optimization calculations, we solved the eigenvalue problem and determined that the first 
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twelve modes were needed to fully cover the frequency range of interest. 

Simply Supported Rectangular Plate . — The first optimization of the plate was done for simply supported 
edges and excitation by a central unit harmonic force. The exciting frequency was 370 Hz, which is slightly lower 
than the natural frequency of mode (3,1), 376 Hz. 

Before executing the optimization, the frequency response of the plate with initial constant thickness was 
evaluated at every 10 Hz up to 1000 Hz considering only the first twelve modes. The result is shown by the broken 
line in figure 30 with the abscissa indicating the excitation frequency. The peak near 380 Hz is the resonant response 
of mode (3,1), and therefore, the energy at the excitation frequency of 370 Hz is high, about 0.338xl0 -4 J. Next, the 
optimized shape was calculated. The convergence of the vibration energy as a function of iteration number is shown 
in figure 31, with the response energy shown as a ratio to that of the first iteration. The energy decreases sharply on 
the first iteration and converges to about l/1000th of the initial energy after eight iterations. The frequency response 
of the optimized design is represented by the solid line in figure 30. The location of mode (3,1) shifts to about 
520 Hz, thus greatly reducing the vibration energy at the exciting frequency. 

The optimum shape of the plate is illustrated in figure 32. The major change is an increase in thickness of 
the elements along the longer center line. The area seems to function as a stiffener, which leads to the reduction of 
modal vibration (3,1). The decrease in thickness at the nodal line that was observed in the optimization of the simply 
supported beam (fig. 16) is not observed in this simply supported plate example. 

A second optimum shape was calculated, in this case for an exciting force at 390 Hz, which is a little higher 
than the natural frequency of mode (3,1). The calculated optimal shape is illustrated in figure 33. The frequency 
response and the convergence of energy are shown in figures 34 and 35, respectively. In this example, the thickness 
increases at the central portion of the plate rather than along a line as in the previous example. The central portion 
appears to function as a mass, and the natural frequency of mode (3,1) is shifted by the optimization to a lower fre- 
quency, about 314 Hz (fig. 34). 

Comparing the frequency responses of the two optimized shapes (figs. 30 and 34), the energy at the excita- 
tion frequency is higher for the case of f ex = 390 Hz. This demonstrates that the minimized energy obtained for this 
excitation is a local minimum. These results also indicate that in this example, for the purpose of reducing vibrations 
in a frequency band that includes the mode under consideration, a shift of a natural mode of vibration to a higher- 
frequency region is preferred, rather than a shift to a lower-frequency region. 

Clamped Rectangular Plate . — As a second example of optimization of plates, the same geometry as the 
previous example was used except the edges were modeled as clamped rather than simply supported. As in the previ- 
ous example, the plate was excited by a central unit harmonic force. The exciting frequency used was 530 Hz, which 
is a little lower than the natural frequency of mode (3,1) at 542.6 Hz. The frequency response of the plate with initial 
thickness of 5 mm is shown by the broken line in figure 36. The energy at 530 Hz was then minimized by the opti- 
mization procedure. The convergence of the energy is shown in figure 37. The energy decreases rapidly on the first 
iteration and is considerably close to the converged value. The frequency response of the optimized plate is shown 
by the solid line in figure 36 and the peak of mode (3,1) shifts to about 680 Hz. 

The optimum shape of the plate is illustrated in figure 38. The thicknesses of elements along the longer 
center line were increased considerably, acting as a stiffener to control the vibration of mode (3,1). However, the 
shape is not as conspicuous as compared to the increase in thickness for the optimum shape of the simply supported 
plate. The thickness of the optimized clamped plate is small at the edges, similar to the effect seen for the clamped 
beam solution. 

The clamped plate was again optimized, this time for an exciting force with a frequency of 550 Hz. The 
optimum shape obtained is shown in figure 39. Since the exciting frequency is slightly higher than the natural fre- 
quency, the optimization shifts the natural frequency to a lower-frequency region in the same way as the simply 
supported plate. The thickness increases at the central portion of the plate and appears to have the function of a tun- 
ing mass. This characteristic is more clear in comparison with the optimum shape in figure 33. Figures 40 and 41 
show the frequency response and convergence of energy, respectively, for this second optimization of the clamped 
plate. The optimization needs more steps to obtain the solution compared to the case of simply supported edges. 

In comparing the results for this example shown in figure 40 with those of figure 36, we conclude that the 
shift of natural frequency to a lower-frequency region is slightly better than the shift to a higher-frequency region if 
the objective is just limited to minimize the energy near the natural frequency of 542.6 Hz. This is opposite to the 
result discussed in the simply supported plate section. However, the difference of the energies is small in this case. 
Therefore, other considerations might be considered in this case for selecting one optimized design over the other. 
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In this chapter, the developed computer code for optimization was applied to the design of beams and 
plates. Clear reductions in vibration through optimization was confirmed by these studies. The optimal shapes have 
areas of increased thickness that act as either a stiffener or a mass. If the frequency of vibration to be minimized is 
slightly lower than the frequency of a nearby natural mode, then this method shifts the frequency of that mode to a 
higher-frequency region. Consequently, since the opposite is also true, this fact can be used to control the direction 
of the natural frequency shift. 

The optimized design should be examined as to whether it is the local minimum. In some cases, shifting the 
natural frequency in the opposite direction can yield further significant reductions in vibration. The results of this 
method are somewhat similar to the usual method of avoiding resonance by shifting the natural frequency through 
strategically locating stiffeners or masses. However, the advantage of using this optimization code is that the user 
need not design or locate the stiffener or mass. The procedure both calculates the shifting of the natural frequency 
and also automatically designs and locates the stiffener or mass. 


REDUCTION OF GEARBOX VIBRATION 
Modeling and Eigenvalue Analysis of Gearbox 

The optimization procedure was applied to the experimental gearbox of the gear noise rig at NASA Lewis 
(ref. 5). The origin of a global X, Y, Z coordinate system was defined at the center of the bottom plate with the Y axis 
parallel to the gear shafts and with the Z axis normal to the bottom plate. The wall center distances of the gearbox 
were approximately X = 324 mm, Y = 248 mm, and Z = 270 mm. The box was modeled using 256 triangular shell 
elements. Although more elements would be desired for a more precise analysis, the number of elements was fixed 
at 256 to limit the required computing resources. 

The pattern of the finite-element mesh is shown in figure 42. This gearbox design includes stiffeners and 
four circular holes for bearings. These structural details, however, were neglected in this model. The mesh is sym- 
metrical about the X axis but not symmetrical about the Y and Z axes. This was done so that the bearing center loca- 
tions would coincide with a nodal point. For boundary conditions, the four bottom corners of the gearbox were fixed 
in all directions. The material constants used were modulus of elasticity E = 206 GPa, Poisson's ratio v = 0.3, mass 

o 1 

density = 8000 kg/m , and coefficients for defining the proportional damping matrix a = 1 .0 s and 
P = 5.0 x10“ 7 s. 

The six surfaces were identified by the numbers 1 through 6 as follows: surface l is the top, 2 is the plate at 
Y = 124, 3 is the plate at X = 162, 4 is the plate at Y = -124, 5 is the plate at X = -162, and 6 is the bottom. Sur- 
faces 2 and 4 support the shafts. The thicknesses of the elements near the bearing locations (—121 < X < 82 and 66.5 
< Z < 193.5) were 22.4 mm with the other elements being 6.3-mm thick. 

Before calculating the optimum design, the eigenvalue problem was solved to determine the dynamic char- 
acteristic of the gearbox (see table IV for the result of the eigenvalue study). The mode shape is represented as (m x , 
m y ) indicating that (m+1) and (m+1) nodal lines are observed in X and Y directions. Sketches of the mode shapes are 
shown in figure 57. Comparing the calculated results to experimental measurements (ref. 1), the calculated funda- 
mental natural frequency at 475 Hz is about 4.4 percent less than the measured value. Though other measurement 
data is available in reference 1, none of the other measured modes coincide with the results of the calculations. This 
may be a result of neglecting the effects of stiffeners and shafts in the structural model used for this research. 


Vibration of Gearbox Caused by Harmonic Exciting Forces 

Unit-exciting forces were used in this research. Since the gear system’s exciting frequencies may be esti- 
mated by spectral analysis of the force, the problem may be considered as harmonic forced vibration of the gearbox. 

The gearbox vibration response to harmonic excitation was computed for unit harmonic forces applied at 
the bearing positions in the direction of the gearing line of action. The total vibration energy was calculated at every 
10 Hz in the region of 500 to 1500 Hz. In addition, the contribution of each gearbox surface to the total gearbox 
energy was computed over the same frequency range (see fig. 43). Ten peaks of resonance are seen in the figure at 
approximately 540, 580, 610, 770, 940, 1240, 1280, 1320, 1400, and 1460 Hz. The figure also shows the magnitude 
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of vibration energy at these resonance frequencies and the contribution of each surface. For example, it is clear that 
the vibration energy at 940 Hz is mainly caused by the vibrations of the top and bottom plates. 

Of course, the mode shape obtained by the eigenvalue analysis gives only the shape and not the magnitude 
of the vibration response. In other words, it represents the relative displacement at the nodes. However, all the eigen- 
vectors calculated in this analysis were normalized with respect to the mass matrix. Therefore, if the component 
related to a certain node is picked up from every mode shape, it might represent the contribution of the mode to the 
relative magnitude of displacement. The components of displacement at the positions of the bearings were collected 
from the results of the calculations and normalized with respect to the maximum value to obtain relative displace- 
ments. The results, normalized as the ratio to the maximum value, are illustrated in figure 44. For example, figure 44 
(a) shows both the displacement in the direction of the X axis and the rotation about the Z axis at the position of 
bearing 1 . The abscissa indicates the order of mode from 1 to 30. The figure shows that the maximum displacement 
in the direction of X axis originates from the 27th mode shape. 

The displacements uj, wj, U 2 , and W 2 were defined as u and w being the displacements in the direction of 
X and Z axes with subscripts 1 and 2 indicating the positions of bearings 1 and 2, respectively. The resonant frequen- 
cies were then obtained from figure 44. In case of spur gears, the work done, U, by the unit load, which is applied at 
the position of bearings along the line of action, can be estimated by the following expression 

U = |uj — u-, sin a + |wj - W 2 cos a (55) 

where a is the operating pressure angle. The work done is also the energy input into the gearbox by the exciting 
force. Therefore, it should have a strong correlation to the vibration energy of the gearbox. Figure 45 shows this 
comparison for the ten resonant frequencies. The vibration energy is approximately proportional to the calculated 
work done even though the calculation is not the actual value but the relative work. This figure suggests that the 
magnitude of gearbox vibration is sufficiently estimated from the relative work, i.e. from the results of the eignvalue 
analysis. 


Optimum Design of Gearbox 

Top Plate of Gearbox . — The gearbox considered here is made of four side plates and a bottom plate welded 
together along with a top plate that is bolted in place. Therefore, strictly speaking, the gearbox model should con- 
sider and include the function of the bolted joints. However, such a model would have a very complex formulation. 
To simplify the formulation, the bolted connections were neglected, and the gearbox was modeled as a simple box 
with all surfaces continuously connected. 

First, the eigenvalue problem for the gearbox was analyzed, and the first 30 modes were selected for use in 
the forced-vibration analysis. Next, the forced-vibration response for the gearbox with all plates at their initial thick- 
ness was calculated. The response was calculated in the range of 400 to 2400 Hz for unit harmonic exciting forces 
applied at the bearing locations and in the direction of the gearing line-of-action. The calculated response in terms of 
vibration energy is shown by the broken line in figure 46(a), and details of the response in a 500- to 
1500-Hz range are shown in figure 46(b). Next, the optimal shape of the top plate was calculated to minimize the 
response to an exciting frequency of 1260 Hz. This frequency coincides with the gear mesh frequency when the test 
gears of 28 teeth operate at 2700 rpm. The vibration energy of the initial design at this frequency is about 
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0.976x10 J. This energy decreases to about 1/5 of this initial value after six iterations of optimization as shown in 
figure 47. This solution was executed on a CRAY XMP/28 computer system, taking approximately 34 sec of cpu to 
compile with the CFT77 compiler and 875 sec of cpu for completing the optimization. The frequency response after 
optimization is shown by the solid line in figure 46. The peak of resonance at 1280-Hz moves to about 1320 Hz. 

Assessing the entire frequency range, the variation of the response looks small. However, the energy is 
greatly reduced at 1200 to 1300 Hz. The optimum shape of the top plate is illustrated in figure 48. The thicknesses of 
the elements increase above bearing 1 and near the edges that connect to surfaces 2 and 4. From figure 44(b), the 
displacement along the Y axis is large in the 12th mode at 1277 Hz. The increase in thickness tends to control this 
displacement by increasing the rigidity of rotation around the X axis. 
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Figure 49 shows the energy share of each surface before and after the optimization. The energy reduction of 
the top plate is relatively small, however, the changes in the top plate design also produces vibration reductions of 
the other surfaces. The frequency response of the energy of each surface is shown in figure 50. In these figures, the 
energy of each surface is compared with the overall energy of the gearbox after optimization. It is clear from the 
results that the peak of resonance at about 1320 Hz is caused by the vibration of surfaces 2 and 4, and that the new 
resonance at 1090 Hz is mainly caused by the vibration of the top plate. 

Design of Top Plate and Two Side Plates of Gearbox . — Next, the two side plates and the top plate were 
optimized simultaneously. Side plates 3 and 5 were selected to be varied by optimization, as they do not support the 
shafts. The condition of excitation used was the same as the previous example, and again, the first 30 modes were 
considered in the analysis. 

The frequency response of the gearbox with initial thickness is shown by the broken line in figure 5 1 . The 
overall vibration energy of the box was minimized for an exciting frequency of 1260 Hz by allowing the element 
thicknesses of the top and two side plates to vary. The frequency response after optimization is shown by the solid 
line in figure 51. The energy is reduced considerably in the region of 1 180 to 1370 Hz. The reduction of vibration 
energy is larger than the first example, where only the design of the top plate was varied. 

The optimization converges such that the response energy is about 1/30 of the initial value as shown in 
figure 52. The optimum shapes of the top plate and two side plates are illustrated in figure 53. The thickness distri- 
bution of the top plate is similar to the top plate of the first example; however, the weight of the top plate is reduced 
about 15 percent by the optimization as compared to the initial weight. The optimum shapes of plates 3 and 5 are 
fairly close to each other. The thicknesses increase toward the bottom corners of the side plates. The weight of plate 
5 is reduced by 7 percent while the weight of plate 3 increases about 23 percent. The total weight of the gearbox was 
constrained to be equal. 

Figure 54 shows the energy share of each surface before and after optimization. The energy reduces for 
every surface, but the reduction for plate 3 is most remarkable. The frequency response of the energy in each surface 
is shown in figures 55 and 56. Figure 55 shows the contribution of each surface on the total energy, and figure 56 
demonstrates the change of frequency response of each surface by the optimization. 


CONCLUDING SUMMARY 

A method was developed to reduce gearbox noise via structural optimization of the housing for low vibra- 
tion. The method is based on the theories of finite elements, modal analysis, and optimization techniques. A three- 
node, triangular- shell finite element with 18 degrees of freedom was developed and implemented for structural and 
modal analysis. 

For optimization, the element thickness was the design variable, while the overall vibration energy in 
response to unit harmonic excitation was the objective function to be minimized. The vibration-response energy 
sensitivity with respect to the design variable was derived and expressed by the sensitivities of eigenvalues and 
eigenvectors. The optimal designs were constrained to have a constant total weight for the structure. Optimal solu- 
tions were computed by the gradient projection method with a unidimensional search procedure. The computer pro- 
gram, REVISE, was written to perform the calculations. 

The finite-element and modal analysis parts of REVISE were verified by comparing the numerical solu- 
tions to exact analytical solutions for (1) the static analysis of a square plate, cantilever beam, and cylindrical shell, 
and (2) the eigenvalue analysis of a rectangular plate. Good agreement between the numerical and analytical solu- 
tions were obtained. 

The optimization procedure was demonstrated by calculating the optimum shapes for minimum vibration 
response for simply supported and clamped beams and plates. The optimal shapes had areas of increased thickness 
that acted as either a stiffener or mass. Forced-response analysis of the optimized shapes clearly demonstrated that 
significant vibration reductions were achieved via the optimization procedure. It was demonstrated that if the vibra- 
tion frequency to be minimized was slightly lower than a nearby mode, the optimizer will shift the frequency of that 
mode to a higher region. The opposite is also true. It was also shown that optimized designs should be examined to 
determine whether the minimum obtained is a local one since, in some cases, further significant reductions in vibra- 
tion energy can be obtained by shifting the frequency in the opposite direction. The optimizer produces similar ef- 
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fects to those of the usual procedure of shifting resonant frequencies by strategically adding a stiffener or mass to the 
design. However, an advantage of using the optimizer is that the stiffener or mass is designed and located 
automatically. 

The computer code was applied to optimize the housing of the NASA Lewis Research Center’s gear noise 
rig test gearbox. Two solutions were obtained: once changing the shape of the top plate only and a second time chan- 
ging the shape of the top and two side plates simultaneously. The shapes were optimized to minimize the vibration 
energy in response to unit harmonic forces applied at the bearing locations along the gearing line of action. 

The total weight of the structure was constrained to be constant. The solutions of the forced-vibration prob- 
lem of the gearbox showed that the vibration energy was approximately proportional to the relative work, which was 
calculated from the bearing forces and the relative displacements estimated from the mode shapes. This suggests that 
the magnitude of forced vibration is approximately estimated from the results of eigenvalue analysis. Compared to 
the initial design, optimizing the top plate only reduced the housing vibration energy by about 1/5. The contribution 
of vibration of each surface to the total vibration was calculated. By viewing the results this way, it was clear that 
optimizing the top plate not only reduced the vibration of that surface but also of the other surfaces as well. 
Futhermore, simultaneously optimizing the shape of the top and two side plates reduced the housing vibration en- 
ergy by about 1/30 compared to the initial design. 
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TABLE I.— DEFLECTION OF SIMPLY 
SUPPORTED SQUARE PLATE 


Mesh 

dimension 

Loading condition 3 

Concentrated, 

W m ax/(PL 2 /D) 

Distributed, 

w m ax/(PL 4 /D) 

Finite element mesh 

8x8x2 

12x12x2 

16x16x2 

0.011636 

.011635 

.011633 

0.004087 

.004081 

.004076 

Closed-form solution 


0.01160 

0.004062 


a Flexural rigidity of plate denoted by D. 


TABLE II.— NATURAL FREQUENCIES OF SIMPLY SUPPORTED 
RECTANGULAR PLATE 


Number of 
nodal lines, 
n 

Mesh 

dimension 

Frequency, Hz 

Number of 
nodal lines, 
m 

i 

2 

3 

4 

5 

6 

i 

4x4x2 

103.0 

198.0 

380.3 

659.8 




6x6x2 

104.9 

204.1 

376.1 

643.2 

1007 



8x8x2 

105.7 

205.0 

376.1 

624.2 

969.1 

1414 


analy. 3 

108.9 

209.4 

376.9 

611.4 

912.9 

1281.4 

2 

4x4x2 

310.3 

383.2 

539.2 





6x6x2 

317.3 

404.0 

558.2 

804.9 

1130 



8x8x2 

321.2 

412.4 

573.4 

806.6 

1141 

1567 


analy. 3 

335.0 

435.5 

603.0 

837.5 

1139.1 

1507.6 

3 

6x6x2 

669.6 

742.6 

864.0 

1099 

1391 



8x8x2 

677.4 

756.9 

900.8 

1115 

1426 

1824 


analy. 3 

711.9 

812.4 

979.9 

1214.4 

1515.9 

1884.5 

4 

6x6x2 

1169 

1228 

1359 

1555 




8x8x2 

1173 

1243 

1370 

1538 

1852 

2216 


analy. 3 

1239.6 

1340.1 

1507.6 

1742.1 

2043.6 

2412.1 

5 

8x8x2 

1810 

1874 

1992 

2180 




analy. 3 

1918.0 

2018.5 

2186.0 

2420.5 

2722.0 

3090.5 


a Closed-form analytical solution. 


TABLE III.— VARIATION OF ELEMENT THICKNESS OF SIMPLY 
SUPPORTED BEAM 


Element 





Element thickness, 

mm 




Iteration 

0 

1 

2 

3 

4 

6 

8 

10 

12 

14 

15 

i 

5.00 

3.63 

2.95 

2.68 

1.96 

1.42 

1.00 

1.00 

1.00 

1.00 

1.00 

2 

5.00 

3.66 

2.98 

2.72 

2.65 

4.24 

4.36 

4.46 

4.49 

5.72 

5.56 

5 

5.00 

4.18 

4.53 

4.59 

5.26 

6.83 

7.16 

7.39 

7.58 

9.28 

9.44 

6 

5.00 

4.11 

4.41 

4.44 

3.47 

1.95 

1.46 

1.00 

1.00 

1.00 

1.00 

9 

5.00 

4.93 

5.36 

5.37 

3.74 

1.04 

1.00 

1.00 

1.00 

1.00 

1.00 

10 

5.00 

5.01 

5.49 

5.56 

6.89 

8.52 

8.84 

9.07 

9.28 

10.00 

10.00 

13 

5.00 

5.86 

5.97 

6.09 

7.65 

9.26 

9.70 

10.00 

10.00 

10.00 

10.00 

14 

5.00 

5.79 

5.86 

5.93 

4.95 

2.58 

1.93 

1.35 

1.00 

1.00 

1.00 

17 

5.00 

6.40 

6.21 

6.29 

6.22 

5.69 

5.53 

5.31 

4.65 

1.00 

1.00 

18 

5.00 

6.43 

6.24 

6.34 

7.20 

8.48 

9.00 

9.42 

10.00 

10.00 

10.00 
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TABLE IV.— NATIONAL FREQUENCIES AND MODE SHAPES OF 
NASA TEST GEARBOX 


Order 

of 

eigenvalue 

Frequency 

ft), 

Hz 

Surface number, nodal lines 

l(m x , m y ) 

2(m x , m z ) 

3(m y , m z ) 

4(m x , m z ) 

5(m y , m z ) 

6(m x , m y ) 

i 

475 

+ (U) 

“(1,1) 

+ (1,1) 

(U) 

+ (U) 

+ (1,1) 

2 

540 

(1,2) 

11,1) 

(2,1) 

+ (U) 

(2,1) 

0,2) 

3 

543 

“(1,1) 

0,2) 

(1,2) 

(1,2) 

0,2) 

+ (1,1) 

4 

582 

+ (1,1) 

+ G,1) 

-(U) 

+ (U) 

-(U) 

+ (U) 

5 

607 

(2,1) 

(2,1) 

+ (U) 

(2,1) 

-(U) 

(2,1) 

6 

768 

+ (1,1) 

+ (U) 

+ (U) 

+ (1,D 

+ (U) 

+ (1,1) 

7 

942 

(2,1) 

°(1,2) 

(1,2) 

°(1>2) 

0,2) 

(2,1) 

8 

994 

(2,1) 

(2,1) 

+ (1,1) 

(2,1) 

-(1,1) 

(2,1) 

9 

1041 

0,2) 

°(1,2) 

(2,1) 

°(1>2) 

(2,1) 

0,2) 

10 

1101 

(1,2) 

°(1,2) 

(2,1) 

°(1>2) 

(2,1) 

0,2) 

11 

1239 

(3,1) 

(1,2) 

(1,2) 

(1,2) 

0,2) 

(3,1) 

12 

1277 

(2,1) 

(2,1) 

(1,2) 

(2,1) 

d (U) 

d (2,l) 

13 

1321 

(2,1) 

(2,1) 

(1,2) 

(2,1) 

0,2) 

(2,1) 

14 

1326 

0,2) 

d ( 1 ,2) 

(2,1) 

d (l,2) 

(2,1) 

0,2) 

15 

1401 

c (2 ) 

0,2) 

(1,2) 

(1,2) 

0,2) 

c (2 ) 

16 

1406 

(1,2) 

(3,1) 

(2,1) 

(3,1) 

(2,1) 

0,2) 

17 

1464 

(2,1) 

°(3,1) 

(1,2) 

°(3,1) 

0,2) 

(2,1) 

18 

1520 

(2,2) 

(2,1) 

(2,1) 

(2,1) 

(2,1) 

(2,2) 

19 

1538 

0,2) 

d (2,l) 

(2,1) 

d (2,l) 

(2,1) 

0,2) 

20 

1646 

(2,2) 

(1,2) 

(2,2) 

0,2) 

(2,2) 

(2,2) 

21 

1699 

(3,1) 

d (3,l) 

(1,3) 

d (3,l) 

0,3) 

(3,1) 

22 

1788 

(2,2) 

d (l,2) 

(2,2) 

d (l,2) 

(2,2) 

(2,2) 

23 

1806 

(3,1) 

°(3,1) 

(1,2) 

°(3,1) 

0,2) 

(3,1) 

24 

1866 

(2,2) 

°(2,1) 

(2,1) 

°(2,1) 

(2,1) 

(2,2) 

25 

1923 

(2,2) 

°(1,3) 

(2,2) 

°G,3) 

(2,2) 

(2,2) 

26 

1953 

(3,2) 

°(1,3) 

(2,2) 

°G,3) 

(2,2) 

(2,2) 

27 

2084 

(4,2) 

d (2,l) 

(2,2) 

d (2,l) 

(2,2) 

(2,2) 

28 

2112 

d ( ) 

°(1,2) 

(1,2) 

°(1,2) 

0,2) 

d ( ) 

29 

2347 

d (4,l) 

(3,1) 

(1,3) 

(3,1) 

0,3) 

d (4,l) 

30 

2448 

(3,2) 

°(1,3) 

(2,2) 

°(1>3) 

(2,2) 

(3,2) 


c Circular mode shape. 
d Deformed mode shape. 
°Oblique nodal line. 
+ Convex outward. 
“Concave outward. 
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Node 3 (x 3 , y 3 ) 



Figure 1 . — Shell element in global coordinate system. 


Function 


Data input 


— - RDATA 


Analyzer 

— ELMMAT 

— LOCAL 

— MASMAT 

— STFMAT 

(transform to local coordinate) 
(form element mass matrix) 
(form element stiffness matrix) 


— STATIC 

— STFBND 

— TRANS 

— LOADNG 

— BCBAND 

— CHOLSK 

— STRESS 

(form global stiffness matrix) 
(transform to global coordinate) 
(form the load vector) 

(introduce boundary conditions) 
(calculate displacements) 
(calculate stresses) 


— OBJ FUN 

— STSSQR 

— EIGEN 

— DAMPNG 

— ENERGY 

(form global mass and stiffness) 

— CHOL (calculate eigenvalues and eigen 

— HOBS vectors up to given order) 
(calculate modal damping ratio) 

(calculate vibration energy) 

Optimizer 

— SENSTV 

— STFMAS 

— SENEIG 
— SENDMP 
— SENENG 

(prepare mass and stiffness) 
(calculate sensitivity of eigenvector) 
(calculate sensitivity of damping) 
(calculate sensitivity of energy) 


GPMACC (introduce active constraints) 

— GPMDIR (calculate feasible direction) 

— SEARCH (search the optimum values) 

Output WRHEAD 

— WRSTRC 
— WRSTAT 
— WREIGN 
— WROBJF 
— WRRESP 
— WRSENS 

Figure 2. — Structure of computer program REVISE. Capital names represent sub- 
routine programs. 
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E = 3x1 0 6 psi 
v = 0 

Dead weight 




2x1 0 3 



<p, degree 

Figure 3. — Deflection and circumferential bending 
moment of cylindrical shell by finite element and 
analytical solutions. 



Figure 4. — Error of eigenvalue analysis of simply supported 
rectangular plate for varying mesh density. 



Figure 5. — Pattern of mesh used for optimum design of 
beam. Length L = 200 mm, width B = 10 mm, initial 
thickness tg = 5 mm. 
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Figure 6. — Frequency response of simply supported 
beam with constant thickness across width. Excitation 
frequency f = 270 Hz. 



Figure 7. — Convergence of vibration energy by iteration 
number for optimization of simply supported beam with 
constant thickness across width. Excitation frequency 
f ex = 270 Hz. 



Figure 8. — Optimum shape of simply supported beam with 
constant thickness across width. Excitation frequency 
f ex = ^Op- 



position, x/(L/2) 

Figure 9. — Comparison of calculated optimum shape for 
minimum vibration energy and shape of uniform strength 
beam. 



Figure 10. — Optimum shape of simply supported beam with 
varying thickness across width. Excitation frequency 
f e x = 270 Hz. 
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Ratio of vibration energy Response, J/N 



Figure 11. — Frequency response of simply supported 
beam with varying thickness across width. Excitation 
frequency f = 270 Hz. 



Figure 12. — Convergence of vibration energy by iteration 
number for optimization of simply supported beam with 
varying thickness across width. Excitation frequency 
f ex = 270 Hz. 


0.6x1 0 -4 



Figure 13. — Displacements due to harmonic excitation 
(initial; uniform thickness, optimum 1, beam in figure 8, 
optimum 2, beam in figure 10). 



Frequency, Hz 

Figure 14. — Frequency response of simply supported 
beam with constant thickness across width. Excitation 
frequency f = 2680 Hz. 
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0123456789 10 

Number of iterations 


Figure 15. — Convergence of vibration energy by iter- 
ation number for optimization of simply supported 
beam with constant thickness across width. Exci- 
tation frequency f = 2680 Hz. 



Figure 16. — Optimum shape of simply supported beam with 
constant thickness across width. Excitation frequency 
f ex = 2680 Hz. 



Figure 17. — Optimum shape of simply supported beam with 
varying thickness across width. Excitation frequency 
f ex = 2680 Hz. 



0 1000 2000 3000 4000 5000 

Frequency, Hz 


Figure 18. — Frequency response of simply supported 
beam with varying thickness across width. Excitation 
frequency f = 2680 Hz. 



0 2 4 6 8 10 12 14 16 

Number of iterations 


Figure 19. — Convergence of vibration energy by iter- 
ation number for optimization of simply supported 
beam with varying thickness across width. Excitation 
frequency f = 2680 Hz. 
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Displacement, 


0.5x1 0“ 5 
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0.0 



Initial 

Optimum 1 
Optimum 2 



' 0.0 0.2 0.4 0.6 0.8 

Position, x/( L/2) 


1.0 


Figure 20. — Displacements due to harmonic excitation 
(initial; uniform thickness, optimum 1, beam in figure 
16, optimum 2, beam in figure 17). 



Frequency, Hz 


Figure 21 . — Frequency response of clamped beam with 
constant thickness across width. Excitation 
frequency f = 650 Hz. 



Figure 22. — Convergence of vibration energy by iter- 
ation number for optimization of clamped beam with 
constant thickness across width. Excitation frequency 

f = 650 Hz. 
ex 



Figure 23. — Optimum shape of clamped beam with constant 
thickness across width. Excitation frequency, f = 650 Hz. 
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Initial 

Optimized 



Position, x/(L/2) 

Figure 24. — Comparison of calculated optimum shape 
for minimum vibration energy with shape of uniform 
strength beam. 



Figure 25. — Optimum shape of clamped beam with varying 
thickness across the width. Excitation frequency 
f ex = 650 Hz. 



0 500 1000 1500 2000 


Frequency, Hz 

Figure 26. — Frequency response of clamped beam with 
varying thickness across width. Excitation frequency 

f ex = 650Hz ' 



Figure 27. — Convergence of vibration energy by iter- 
ation number for optimization of clamped beam 
with varying thickness across width. Excitation fre- 
quency f = 650 Hz. 
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Figure 28. — Displacements due to harmonic excitation 
(initial; uniform thickness, optimum 1, beam in figure 
23, optimum 2, beam in figure 25). 



a 


Figure 29. — Pattern of finite element mesh used for optimum 
plate design (a = 600 mm, b = 400 mm, initial thickness 
t = 5 mm). 



Figure 30. — Frequency response of simply supported plate. 
Excitation frequency f = 370 Hz. 
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Figure 31 . — Convergence of vibration energy by 
iteration number for optimization of simply sup- 
ported plate. Excitation frequency f ex = 370 Hz. 


Optimized 



Figure 32. — Optimum shape of simply supported plate. 
Excitation frequency f = 370 Hz. 



Figure 34. — Frequency response of simply supported plate. 
Excitation frequency f = 390 Hz. 



Figure 33. — Optimum shape of simply supported plate. 
Excitation frequency f = 390 Hz. 
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Figure 35. — Convergence of vibration energy by iteration 
number for optimization of simply supported plate. Exci- 
tation frequency f ex = 390 Hz. 



Figure 36. — Frequency response of clamped plate. Exci- 
tation frequency f = 530 Hz. 



Figure 38. — Optimum shape of clamped plate. Excitation 
frequency f ex = 530 Hz. 




Figure 39. — Optimum shape of clamped plate. Excitation 
frequency f ex = 550 Hz. 


Figure 37. — Convergence of vibration energy by 
iteration number for optimization of clamped 
plate. Excitation frequency f ex = 530 Hz. 
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Ratio of vibration energy Response, J/N 



Figure 40. — Frequency response of clamped plate. 
Excitation frequency f = 550 Hz. 



Figure 41 . — Convergence of vibration energy by iteration 
number for optimization of clamped plate. Excitation 
frequency f ex = 550 Hz. 



Y = 124 



ii 


X 

Figure 42. — Finite element mesh pattern used for the 
optimum design of NASA test gearbox. 
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Response, J/N Response, J/N 



500 1000 1500 500 1000 1500 

Frequency, Hz Frequency, Hz 


Figure 43. — Frequency response of NASA gearbox with initial thickness showing total vibration energy response and contri- 
bution of each plate to total, (a) Contribution of plate 1 (top) on vibration energy, (b) Contribution of plate 2 on vibration energy, 
(c) Contribution of plate 3 on vibration energy, (d) Contribution of plate 5 on vibration energy, (e) Contribution of plate 6 on 
vibration energy. 
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Response, J/N 


1 



Frequency, Hz 

Figure 43. — Concluded, (e) Contribution of plate 6 on 
vibration energy. 
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Figure 44. — Relative displacements at bearing location as 
function of mode order, (a) Displacement along X axis and 
rotation about Z axis for bearing 1 . (b) Displacement along 
Y axis for bearing 1 . (c) Displacement along Z axis and 
rotation about X axis for bearing 1 . (d) Displacement along 
X axis and rotation about Z axis for bearing 2. (e) Displace- 
ment along Y axis for bearing 2. (f) Displacement along Z 
axis and rotation about X axis for bearing 2. 
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Figure 45. — Relation between vibration energy of gearbox 
and relative work done by exciting force. 




Figure 46. — Frequency response of NASA gearbox with 
shape of top plate optimized for minimum vibration energy. 

(a) Response in frequency range 400 through 2400 Hz. 

(b) Response in frequency range 500 through 1500 Hz. 
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Ratio of vibration energy 



Figure 47. — Convergence of vibration energy by iteration 
number for optimization of NASA test gearbox top plate 
only. 
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Figure 49. — Vibration energy response of each NASA test 
gearbox plate for initial design and after optimizing top- 
plate. 


Figure 48. — Optimum shape of NASA test gearbox top plate 
for minimum vibration energy when optimizing top-plate 
shape only. 
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Response, J/N 



500 1000 1500 

Frequency, Hz 


Figure 50. — Concluded, (e) Plate 6. 



Frequency, Hz Number of iter ation 

Figure 51 . — Frequency response of NASA gearbox with Figure 52 — Convergence of vibration energy by iteration 

shape of top and two side plates optimized for minimum number for optimization of NASA test gearbox top and 

vibration energy. two side plates. 


Top plate 


Side plate 
(plate 5) 


Side plate 
(plate 3) 


Figure 53. — Optimum shape of NASA test gearbox top and two side plates for minimum vibration energy when optimizing 
all three sides simultaneously. 



NASA gearbox plate 


Figure 54. — Vibration energy response of each NASA test 
gearbox plate for initial design and after optimizing top 
(plate 1) and sides (plates 3 and 5) simultaneously. 
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Figure 55. — Frequency response of each NASA gearbox plate with shapes of top (plate 1) and sides (plates 3 and 5) 
optimized simultaneously for minimum vibration energy compared to total response, (a) Plate 1 (top), (b) Plate 2. 
(c) Plate 3. (d) Plate 5. (e) Plate 6. 
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Figure 56. — Frequency response of each plate before and after optimization for case of optimizing top (plate 1) and 
sides (plates 3 and 5) simultaneously for minimum vibration energy, (a) Plate 1 (top), (b) Plate 2. (c) Plate 3 (side), 
(d) Plate 5 (side), (e) Plate 6. 
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Figure 57. — NASA test gearbox computed mode shapes. 
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Figure 57. — Concluded. 
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